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|> ■ ABSTRACT 



Argon abundances have been derived for a sample of B main- sequence stars 
in the Orion association. The abundance calculations are based on NLTE metal 
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line-blanketed model atmospheres calculated with the NLTE code TLUSTY and 
an updated and complete argon model atom. We derive an average argon abun- 
dance for this young population of A(Ar) = 6.66 ± 0.06. While our result is in 
excellent agreement with a recent analysis of the Orion nebula, it is significantly 
higher than the currently recommended solar value which is based on abundance 
measurements in the solar corona. Moreover, the derived argon abundances in 
the Orion B stars agree very well with a measurement from a solar impulsive flare 
during which unmodified solar photospheric material was brought to flare con- 
ditions. We therefore argue that the argon abundances obtained independently 
for both the Orion B stars and the Orion nebula are representative of the disk 
abundance value in the solar neighborhood. The lower coronal abundance may 



1 Based on observations obtained with the Apache Point Observatory 3.5-meter telescope, which is owned 
and operated by the Astrophysical Research Consortium. 
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reflect a depletion related to the FIP effect. We propose a new reference value 
for the abundance of argon in the solar neighborhood, A(Ar) = 6.63 ± 0.10, 
corresponding to Ar/O = 0.009 ± 0.002. 

Subject headings: astrochemistry — stars: abundances, early-type — Sun: abun- 
dances 



1. INTRODUCTION 



Establishing the abundance of noble gases in the Sun has been a notoriously complex 
task, with significant uncertainties still remaining in some cases. These difficulties primarily 
ensue from the atomic structure of noble gases, in particular their high first ionization po- 
tential and the large excitation energies of the excited levels in neutral atoms. Noble gases 
are therefore in the neutral ground state to a great extent in the solar photosphere and, 
consequently, no spectral lines of these elements have been identified in the photospheric 
solar spectrum. Well-established model stellar atmosphere and quantitative spectroscopic 
techniques therefore cannot be applied to study noble gases in the Sun. Furthermore, these 
species are volatile. Meteoric studies thus yield only a lower limit of the actual abundance 
of noble gases in the solar system and, therefore, are of limited usefulness for inferring the 
solar argon abundance. 

While the solar helium abundance can be determined with a great degree of accuracy 
from helioseismology, abundance studies of heavier noble gases depend on the analysis of 
coronal lines and in-situ measurements of energetic particles from the solar wind. These 
techniques have greater uncertainty than photospheric abundance studies which have been 
applied to many other chemical elements. Moreover, measurements are made relative to a 
reference element, and they actually yield a bundance ratio s such as the Ar/O ratio. In a 
recent compilation of solar abundance works, lAsplund et al.l (120051 ) recommend a solar argon 
abundance, A(Ar) = 6.183 This value significa ntly revises downw ard the e arlier standard 



value (A(Ar) = 6.55. lAnders & Grevessd ll989h : A(Ar) Q = 6.5 0. lLoddersI J2008hh mostly 



because of the down ward revision of the solar oxygen abundance ( lAllende Prieto et al.ll2001 
Asplund et aJl2004h . 



Because of the shortcomings of solar studies, it appears that abundance measurements in 
the solar neighborhood may shed more light on the argon abundance at near-solar metallic- 



3 Abundances are quoted in number densities throughout the paper, in the standard scale where 
A(H) = 12.0. 
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ity. Indeed, this approach is suggest ed by the similar case of neon, for which rec ent analyses 
in B-type stars (ICunha et al.ll2006l ). in the Orion n ebula (jEsteban et al.l 12004 ). and in the 
interstellar medium (ISM) toward the Crab nebula (IKaastra et al.l 120071 ). provide very con- 
sistent results, Ne/O = 0.26±0.01, and therefore indicate that the currently adopted solar 
value, (Ne/O)o=0.15, most likely underestimates the actual neon abundance in the solar 
photosphere and in the solar neighborhood. Defining a reliable set of standard chemical 
abundances is an essential undertaking, because abundance studies are a common tool used 
to study physical processes which are revealed by changes in the chemical composition of 
various astrophysical media, for instance including chemical anomalies in the solar corona, in 
the ISM (e.g., depletion into grains, ISM ionization), or in stellar photospheres (e.g., mixing 
of nucleosynthesis products). 

The optical spectra of early-type stars reveal a number of unblended, weak Ar II lines. 
Despite appearing to be the most fitting targets from which to derive the argon abundance in 
the solar neighborhood , very limited abund ance work on argon in early- type stars has been 
carried out to this day. iKeenan et al.l (119901 ) analyzed two Ar II lines in the blue spectrum of 
5 bright B-type stars in the field, assuming Local The rmodynamic Equilibriu m (LTE), and 
they derived an argon abundance, A(Ar)= 6.49 ±0.1. iHolmgren et al.l ( 19901 ) supplemented 
this study with an investigation of the importance of departures from LTE, which were found 
to be small. The NLTE result (A(Ar)= 6.50 ± 0.05) thus does not differ significantly from 
the LTE analysis; these NLTE calculations, however, were based upon non-blanketed model 
atmospheres. 

In this paper we report on the results of a new NLTE analysis of a sample of 10 main- 
sequence B-type stars which are members of the young Orion association. The argon abun- 
dances are derived from NLTE line-blanketed model atmosphere calculations. We argue 
that the argon abundances obtained for this young population may well represent the actual 
argon content of the solar neighborhood and that the Orion B stars can be used as good 
proxies to define the argon abundance in the solar photosphere. Sect. 2 briefly presents the 
new observational material, while the NLTE calculations are detailed in Sect. 3. We dis- 
cuss our results in Sect. 4, comparing them to other determinations of the argon abundance 
in different sites. In our conclusions (§5), we propose a new reference value for the argon 
abundance in the solar neighborhood. 



2. OBSERVATIONS AND DATA REDUCTION 



The sample under study consists of 10 main-sequence early B-type stars, whi ch have 
been previously analyzed in our recent work devoted to neon (ICunha et al.l 120061 ) . New 
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spectra were obtained on the nights of 2007 January 29/30, February 5/6, and March 2/3, 
using the ARC echelle spectrograph on the 3.5 m telescope at the Apache Point Observatory. 
The ARC echelle covers the entire visible spectrum on a single 2048x2048 detector with a 
resolving power of R ~ 35,000. Exposure times for the different targets are listed in Table [TJ 
Because of partly cloudy conditions, some stars were observed repeatedly in order to achieve 
a signal-to-noise ratio of at least 100 in the co-added spectra. 

The data were reduced and the spectra extracted following fairly standard procedures. 
Bias frames and flat field exposures were combined; flat frames included frames with a blue 
filter to insure a high signal-to-noise ratio across the entire spectrum. Individual orders were 
extracted from both the flats and the object frames, after subtraction of a relatively small 
scattered light component. Some care was taken with the extraction parameters because of 
the relatively narrow width of the orders in the spectrograph. The extracted spectra were 
flattened using the one dimensional extracted flats. Wavelength calibration was performed 
with Th-Ar exposures, and applied to the object spectra. Noise spectra were obtained using 
the known gain and readout noise of the detector. The extracted spectra were coadded 
when multiple observations were available. Finally, the individual orders were combined, 
using the noise spectra to properly weight the different orders in regions where there was 
spectral overlap. Sample spectra of all target stars in the region around A 4429 A are shown 
in Fig. [TJ 



3. NON-LTE ABUNDANCE CALCULATIONS 



The NLTE calculations have been performed wit h the NLTE model stellar atmosphere 
code TLUSTY flHubenvl fl988l : iHubenv fc Lan3ll995h fl The TLUSTY model atmospheres 
assume hydrostatic and radiative equilibria, and incorporate all essential sources of opacity 
while allowing for departures from LTE. A large grid of NLTE fully line-blanketed model 
atmospheres of early B-type stars has been made available recently (BSTAR2006; Lanz & 
Hubeny 2007). The present calcu lations follow the BSTAR2006 models, but for larger Ne i-ll 
model atoms (ICunha et al.l 120061 ) and for adding extensive Ar i-iv model atoms (see §3.11) 
in the NLTE calculations. The NLTE argon calculations thus benefit from an extensive and 
detailed NLTE treatment of all background opacity sources. 

We have adopt ed for the new models the same stellar parameters of the target stars as 



Cunha et al.l (120061 ) , see Table [3j The stellar parameters were derived by ICunha fc Lambert 



( 119921 ). who provided a detailed discussion about uncertainties. We examine their impli- 
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cations for the argon abundances in §3.21 Like in BSTAR2006, a constant microturbulent 
velocity of 2 km s _1 was assumed in the NLTE model atmosphere calculations. 

NLTE synthetic spectra were then calculated with the spectrum synthesis code SYN- 
SPEC0 using the NLTE level populations for all levels included in the TLUSTY models. 
We have adopted in this step a microturbulent velocity of 5 km s _1 that is typical for early 
B-type main-sequence stars. We stress, however, that the Ar II lines analyzed in all 10 stars 
are relatively weak and, therefore, the results are quite insensitive to the microturbulence 
parameter (see §3.2j) . 



3.1. Argon model atom 

The argon model atom consists of 71 levels of Ar I, 54 levels of Ar II, 44 levels of 
Ar III, 36 levels of Ar IV, plus the ground state of Ar V. Because the level structure of 
Ar I is better described with intermediate j.K'-coupling (similarly to Ne I), the Ar I model 
atom was constructed to account for fine structure. It includes the lowest 71 fine-structure 
levels (up to the 6s and 6p levels, E < 121,500 cm -1 ). Experimental energy level data 
and li ne oscillator strengt hs were extracted from the Atomic and Spectroscopic Database at 
NIST jMartin et alill999T lPl The photoionization cro ss-sections were taken from the Opacity 
Project (OP) database Topbase ( Cunto et al. 19931 ) Since OP calculations assume LS- 



coupling, we followed a procedure analogous to the case of Ne I (jSeatonlll998l ) for relating 
fine-structure levels to LS'-multiplet levels and thus for assigning the appropriate photoion- 
ization cross-sect i on to each level. The collisional excitation rates were considered using the 
Van Regemorterl (119621) fo rmula, including a modification for neutral atoms (as also used 
by lAuer fc Mihalad (119731 )) ; for col l isiona l ionization the Seaton formula was used (for a 
synopsis of expressions, see iHubenyl (119881 )). 



For the Ar n-iv ions, we closely followed the approach described in lLanz fc Hubeny 



(120031 ) regarding the treatment and inclusion of atomic data in model atmosphere calcula- 
tions. In particular, the bulk of the data are taken from Topbase. The level energies were 
updated with the more accurate experimental values whenever available in the NIST ASD 
database. For a detai led description o f the p rocess of setting up the atomic data, the reader 
should refer to §4 of lLanz fc Hubenyl (120031 ). In particular, we include in each model ions 
all the energy levels which are in Topbase and below the respective ionization limit, that is, 



See http://physics.nist.gov/PhysRefData/ASD/index.html 



See http://vizier.u-strasbg.fr/topbase 
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typically up to n = 10. The Ar II model atom includes individually the lowest 42 levels up 
to Ad 2 D (E < 193, 000 cm" 1 ). The 206 higher levels are grouped into 12 superlevels, 6 in the 
doublet and 6 in the quartet systems. The Ar III model atom includes the lowest 27 levels 
up to 3d 3 P° (E < 215, 000 cm" 1 ). The 349 higher levels are grouped into 17 superlevels, 6 in 
the singlet, 6 in the triplet and 5 in the quintet systems. The Ar IV model atom includes the 
lowest 22 levels up to 4s 2 D (E < 270, 000 cm" 1 ), and the first 3 sextet levels. The 322 higher 
levels are grouped into 10 superlevels, 5 in the doublet and 5 in the quartet systems. Five 
excited sextet levels, below the ionization limit, are grouped into an additional superlevel. 



3.2. Abundances from selected Argon lines 



The optical spectra of early B-type stars contain a number of unblended Ar II lines. 
The sample of 11 Ar II lines analyzed in this st udy is presented in T able [2j The Ar II 
gf- values for these transitions were extracted from iBennett et al.l (119651 ) via the NIST ASD 
database. The accuracy assessment of these oscillator strengths by NIST is very high, with 7 
values estimated to be as accurate as 3% while the 4 other lines have estimated uncertainties 
smaller than 10%. The linelist used in co mputation of t he synthetic spectra around each 
of the Ar II transitions is an update of the iKuruca (I1994J ) linelist and is available from the 
TLUSTY website. 4 

Argon abundances were derived from best fits between synthetic and observed spectra. 
Sample fits for 6 Ar II lines (AA 4426; 4430; 4735; 4764; 4806 and 5062 A) in the star 
HD35299 are shown in Fig. [2j The best-fit synthetic profiles for each line were calculated 
with the argon abundances indicated in the different panels of the figure. The sensitivity of 
the fits to changes in the argon abundance by ± 0.1 dex is illustrated for the Ar IIA4426 line 
in the top left panel of Fig. [2j The final argon abundances for the target stars are given in 
Table [3l They are the averages of the abundances derived from the individual Ar II lines. The 
corresponding standard deviations represent the line-to-line scatter, which may be attributed 
in part to the remaining uncertainties on the /—values (up to 0.05 dex). For completeness, we 
also show in Table [3] the oxygen abundances obtained previously for the target stars, because 
coro nal studies provide the ab und ance ratio Ar/O. T he oxygen abundances are extracted 
from lCunha fc Lambertl (11994 ); see lCunha et al.l (120061 ) for a discussion of consistency checks 
between these oxygen abundances and TLUSTY. Abundances are listed as number densities 
in the standard logarithmic scale, where the hydrogen abundance is A(H)=12. 

The uncertainties in the derived argon abundances can be evaluated by investigat- 
ing abundance changes resulting from model atmospheres computed with modified stellar 
parameters. The error budget in the argon abundances is dominated by the uncertain- 
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ties in T e ff. The adopte d effective temperatures are derived from Stromgren photometry; 



Cunha fc Lambertl (119941 ) carefully assessed the uncertainties to 5T eff = ±3% at most. We ver- 



ified that the adopted T e g- values and uncertainties are consistent with temperature-sensitive 
spectroscopic diagnostics such as the Si IIIAA4553, 4568, 4575 lines. These uncertainties 
translate to an error on the derived abundances, £A(Ar)= ±0.08 dex. On the other hand, 
the uncertainty on gravity, Slog g= ±0.1 dex, has a very limited effect on the abundances, 
#A(Ar)= =f0.01 dex. The errors due to uncertainties in the adopted microturbulent velocities 
are small given that the sample Ar II lines are all reasonably weak. A change in microturbu- 
lence of 1.5 km/s results in <5A(Ar)= ±0.01 dex. The total error from uncertainties in these 
stellar parameters is ±0.08 dex, using a quadratic error summation. 



4. DISCUSSION 



The argon abundance cannot be measured in the photospheres of the Sun and of other 
cool stars. We discuss and compare argon abundances derived in different sites, including the 
solar corona, in-situ measurements in the solar system, the interstellar medium (ISM), H II 
regions and planetary nebulae, along with our analysis of early-type stars. The most repre- 
sentative results are discussed in the next sections. This comparison provides for a reliable 
determination of the abundance of argon in the solar neighborhood, which can then serve as 
a reference point for further st udies. A s i milar overview of argon abundance determinations 
has been recently presented by iLodderd (120081 ). 



4.1. Argon in early-type stars and departures from LTE 

The weighted (as well as straight) average of the NLTE abundances measured in 10 
early B stars of the Orion association, and listed in Tabled is <A(Ar)>= 6.66 ± 0.06. All 
the Ar II lines give a consistent result, as indicated by the line-to- line scatter for each star 
which is always smaller than 0.1 dex, and the abundances span a very tight range (6.59 < 
A(Ar) < 6.76). The error budget from uncertainties in stellar parameters is ±0.08 dex. 
Furthermore, the available /-values are of high accuracy, generally better than 3%. We 
therefore believe that the accuracy of the derived argon abundance is likely of the order of 
0.1 dex or better. 



Keenan et al.l (Il990l ) carried out an early, but very limited, study of argon in 5 nearby 
field B-type stars. They observed two Ar II lines (AA4590, 4658). One star was excluded 
from the sample because only Ar IIA4590 was marginally detected. Argon abundances were 
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derived from equivalent widths and LTE model atmospheres. Results from a second star 
(HR 1350) were also excluded because the derived abundance differed by 0.25 dex relative 
to the 3 other stars, which was attributed to larger uncertainties in the stellar parameters. 
They derived an abundance, A(Ar)= 6.49 ± 0.1, based on 4 measurements only. Their value 
is 0.17 dex lower and barely consistent with our results. However, if we retain the abundance 
derived for HR 1350 that they excluded, and if we correct for small systematic differences in 
the /-values (0.024 and 0.046 dex for AA4590, 4658, respectively), we find a better agreement 
(A(Ar)= 6.54 ± 0.15). The larger scatter seems to be a more likely estimate of the actual 
uncertainties of their results. 



Holmgren et al.l (119901 ) supplemented Keenan et al.'s study and evaluated the impor- 
tance of de partures from LTE. They ap plied the NLTE line formation codes DETAIL and 
SURFACE n Giddinga 1 1 98 ll ; lButlerlll984j ) and found that the corrections from LTE are very 
small. Holmgren et al. derived an argon abundance, A(Ar)= 6.50 ± 0.05. However, their 
NLTE analysis is based on unblanketed NLTE model atmospheres, while Keenan et al. used 
metal line-blanketed LTE models. Their comparison therefore amalgamates two different 
effects, namely direct NLTE effects in argon and differences arising from the changes in 
atmospheric structure between blanketed and unblanketed models. 

To ascertain better the extent of departures from LTE in argon, we have repeated 
our analysis of the target star HD35299 using a Kurucz LTE blanketed model atmosphere 
and LTE spectrum synthesis. We can then make a direct comparison of two consistent 
analyses with and without the assumption of LTE, but based on fully line-blanketed model 
atmospheres in both cases. The LTE abundances obtained for the sample Ar II lines are 
not significantly different from the NLTE values: < 5 (NLTE - LTE)>^ 0.02 ± 0.01 dex, 
supporting Holmgren et al.'s earlier claim that departures from LTE are small. LTE argon 
abundances can also be derived using the same NLTE model atmosphere, but only setting 
argon populations to their LTE values. This second approach may isolate the non-LTE 
effects in argon since the same model structure is adopted. In this case, our results show 
that NLTE effects are also small, the NLTE abundances being larger than the LTE values 
by_~,0.03 dex. These results therefore support the earlier NLTE study of iHolmgren et al. 



(Il990l ). We stress however the importance of performing extensive NLTE calculations when 
aiming at determining abundances with accuracies as high as 0.1 dex or better. Indeed, even 
small corrections to LTE may easily introduce systematic errors of the order of the expected 
accuracy. 

We examine now the reas ons why the departures from LT E are small for Ar II lines while 
quite significant for Ne I lines ( lSigutlll999l ; ICunha et al.ll2006l ). The top panels of Fig. [3] show 
that the typical NLTE ionization shift toward higher ionization is limited. Because of the 
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high ionization potentials of neon and argon ions, the ionization thresholds are found in the 
extreme ultraviolet where B stars have a low flux. In this spectral range, the opacity is 
dominated by hydrogen and helium. The bound-free transition from the ground states of 
neon and argon are optically thick and in detailed radiative balance in the region of formation 
of neon and argon optical lines. Consequently, the contribution of radiative photoionization 
to departures from LTE remains small. Rather than by a global ionization shift, the line 
strengths are then better described by the departure coefficients, 6j = nj/n*, where rii is 
the population of level % and n* its LTE counterpart. The bottom panels of Fig. [3] show 
a very different behavior for lines in the 3s — 3p Ne I and in the 4s — 4p Ar II systems to 
which the analyzed lines belong. In both cases, the optical lines behave mostly like lines in a 
2-level atom because the resonance lines are optically thick at the depth of formation of the 
optical lines. In the neon case, the (red) lines provide a downward channel, overpopulating 
the lower levels, th us resulting in stronger absorption and yielding thus lower abundances 



(ICunha et al.ll2006l ). On the other hand, the Ar + levels are collisionally coupled to the Ar 2+ 
ground state, as indicated by the fact that all Ar + levels show similar departure coefficients. 
In first approximation, the ratio of the line source function to the Planck function behaves 
like S/B oc b u /bi. In the line forming region, our calculations typically give b u /b[ »s 0.98, thus 
yielding only a small correction to LTE. Finally, the Ne I and Ar II resonance lines are formed 
in the outermost atmospheric layers, around tr oss « 10~ 4 , where they provide a downward 
radiative channel, leading to an overpopulation of the ground states and underpopulation of 
the excited states. 



Among hot stars. I Werner et al.l (120071 ) have recently reported the identification of Ar VII A1063 
in very hot central stars of planetary nebulae and (pre-)white dwarfs. Because of large re- 
maining uncertainties on stellar parameters, and because of processes such as nucleosynthesis 
and elemental diffusion occurring in late evolutionary stages, their derived argon abundances 
cannot be used as a reference value. Several stars however show an argon abundance close 
to the value that we derived for B-type stars, hence supporting the idea that there is no 
large production of argon in the previous AGB evolution and that argon remains essentially 
unchanged in these stages. 



4.2. Nebular emission abundances 



Esteban et al.l (120041 ) derived emission line abundances for several elements in the Orion 
nebula. In particular, they obtained A(Ar)= 6.62 ± 0.05 when adopting their preferred 
value for the temperature fluctu ation parameter, t 2 =0.02. The arg on abundance in the 



Orion nebula derived earlier by iPeimbert fc Torres-Peimbertl ( 119771 ) also ranges between 
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A(Ar)=6.6 (with t 2 =0) and A(Ar)=6.7 (for t 2 =0.35). There is therefore a good indication 
of a consensus regarding the nebular argon abundance in Orion. Furthermore, the nebular 
analysis is fully consistent with our results from B stars in the Orion association (A(Ar) = 
6.66 ± 0.06), even with the tight uncertainties of the two studies. 



Argon abundances can also be derived from planetary nebulae (PNe) . IStanghellini et al. 



( 120061 ) present argon abundances for a large sample of galactic PNe and find an average argon 
abundance for the round nebulae (representing the lower mass and probably local population) 
of <A(Ar)> = 6.11. These argon abundances rely on optic a l argo n lines and are deemed 
as more uncertain in their study. iPottasch fc Bernard- Salad (120061 ) give a summary of PN 
abundances from measurements recorded with the Infrared Space Observatory. The Galactic 
chemical gradient is apparent from their dataset when plotting abundances as a function of 
galactocentric distances. Therefore, for comparison with the Orion results, we have retained 
only the PNe which are within 1 kpc from the Sun. The average argon abundance of the 
7 PNe satisfying this criterion is A(Ar) = 6.50l 22- The average value is 0.12 - 0.16 dex 
lower than the results from Orion. However, the scatter is much larger and we cannot 
conclude that there is a s i gnific ant difference. The results of our study therefore support 
Pottasch fc Bernard-Salad (120061 ) conclusion that the argon abundance is essentially the same 
in young objects (B stars, Orion nebula) and in older PNe in the solar neighborhood. 



4.3. ISM absorption studies 



Because of its high ionization potential, argon is mostly in the neutral state in the ISM. 
Argon abundances can be derived from column densities measured from the far-ultraviolet 
resonance lines, Ar I AA1048, 1067. 1 Sofia fc Jenkind (119981 ) detected these Ar I lines towards 
several early-type stars using the interstellar medium absorption profile spectrograph during 
the first ORFEUS-SPAS shuttle mission in 1993 September. They found deficiencies in 
Ar I/H I r anging from 0.2 to 0.6 dex, relative to the Ar /H abundance ratio in the Sun and 
in B stars (lAnders fc Grevessdll989l ; iKeenan et al.lll990l ). They argued that this result does 
not provide an indication of argon depletion onto dust grains, but it rather reflects selective 
argon ionization relative to hydrogen because of the larger Ar I photoionization cross-section. 

With the launch of the Far Ultraviolet Spe c trosco pic Explorer, many more line-of-sights 
in the ISM have been observed. iLehner et al.1 (120031 ) carried out an analysis of 23 line-of- 
sights toward white dwarfs and derived argon abundanc es in the local ISM. Th ey confirmed 
deficiencies of the order of 0.4 dex, as found earlier by ISofia fc Jenkind (I1998I ). Moreover, 
the idea that these deficiencies actually map the local ISM ionization is well supported by 
an increase of the abundance scatter outside the local Bubble. It is therefore essential to 
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have a good reference value of the local argon abundance for ISM ionization studies, but 
Ar I column densities are thus of limited use for accurately determining the reference argon 
abundance in the solar neighborhood. 



4.4. Coronal vs. photospheric abundances 



The standard solar argon abundance recommended by lAsplund et al.l (120051 ). A(Ar) = 
6.18 ± 0.08, has been derived from measured abundance ratios in the solar corona and solar 
energetic particles. However, significant differences in the chemical composition of the solar 
corona relative to the photosphere were noticed very early on. These differences were found 
to correlate with the first ionization potential of the studied elements, the so-called "FIP 
effect". While there is still no agreed theoretical explanation of the FIP effect, the current 
general consensus is that elements with lower FIP (< 10 eV) are enhanced in the corona 
by a factor 4 to 6, while elements w ith higher FIP (> 10 eV) have essentially photospheric 
abundances ( Favata fe Micela 2003h. However, most coronal studies based for instance on 



Ar XIV A4412 (e.g.. lYoung et al. 



1997) or on helium-like argon lines observed in X-ray flares 



(jVeck fc Parkinson! Il98ll ; iPhillips et al.l 120031 ) indicate an argon abundance, typically in a 
range between 6.2 < A(Ar) < 6.4, that is lower than in Orion stars and nebula. The Orion 
results thus support the idea that the solar corona exhibits significant underabundances of 
the two elements with a high FIP (> 15 eV). However, we cannot define the origin of these 
underabundances. The FIP effect might either cover a whole range, from overabundances for 
low FIP elements, to photospheric abundances at intermediate FIP and underabundances 
for elements with the higest FIP, or different physical processes might be in play and the 
behavior of low FIP and high FIP elements might be unrelated. 

Since the neon and argon emission lines observed in the coronal and chromo s pheric 
spectra of the Sun are potentially affected by the FIP effect, iFeldman fc Widingl (120031 ) 
argued that it would be much preferable to observe and analyze unmodified photospheric 
plasma. During a low-altitude impulsive flare, photosph eric plasma raised to flare temper- 
ature was observed by Skylab (IFeldman fc Widingl Il9901 ). From the high measured plasma 
density and a Mg/ O abundance ratio similar to the photospheric value, Feldman & Widing 
established the origin of the emission. They determined the Ar/Mg and Ne/Mg abundance 
ratios, yielding A(Ar ) = 6.57 ± 0.12 and A(Ne) = 8.08 ± 0.10, with A(Mg) = 7.53 from 
Asplund et al.l (120051 ). These values are in excellent agreement with the argon and neon 
abundances measured in B stars and in the Orion nebula. 



Argon abundances can also be measured in th e corona of other co ol stars where the 
FIP effect might not be as prevalent as in the Sun. iMaggio et al.l (120071 ) recently reported 
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on a large study of X-ray bright pre-main sequence (PMS) stars in Orion. They derived an 
average argon abundance of A(Ar) = 6.5 lift??. The larger scatter implies that this result 
remains consistent with the abundance measured in Orion B-type stars and, therefore, does 
not support an underabundance of argon in the corona of Orion PMS stars similar to the 
underabundance found in the solar corona. 



4.5. In-situ measurements in the Solar system 



Argon abundances have also been measured from energetic particles in the solar wind, 
from lunar soils, and meteorites. iReamed (119981 ) contended that the most comprehensive 
measurements of element abundances in the solar corona comes from measurements of ener- 
getic particles that have been ejected during coronal mass ejections . Becaus e of th eir origin, 
we expect however that these abundances to show the FTP effect. IReamed (I1998T) reported 
an abundance ratio Ar/O = 0.0033 ± 0.0002. Usi ng lAsplund et al.l (120051 ) oxygen abun- 
dance, this ratio translates to A(Ar) = 6.18 ± 0.03. Ceruttil ( 19741 ) reported a similarly low 
Ar/O ratio in lunar soils, yielding A(Ar) = 6.26 ± 0.1. Finally, CI ch ondrites reveal strong 
depletion because of the volatile nature of noble gases. For instance, lAsplund et al.l (120051 ) 
list an argon abundance in meteorites that is 6 to 7 orders of magnitude lower than all other 
measurements. 



4.6. Abundances in the Solar neighborhood 

Table H] summarizes the argon abundances measured in different sites in the Solar system 
and in the solar neighborhood, using a broad variety of methods. Excluding meteorites, argon 
abundances span a range of 0.6 dex (a factor of 4), 6.1 < A(Ar) < 6.7. On one hand, solar 
corona and in-situ measurements in the Solar system yield a low value (6.2 < A(Ar) < 6.4, 
or A(Ar) = 6.18 as adopted by Asplund et a!.). On the other hand, Orion B-type stars and 
nebula provide a tight agreement with emerging photospheric material during an impulse 
flare, 6.57 < A(Ar) < 6.66. What is therefore the best value to represent the argon abundance 
in the Solar system? Can we use Orion stars and nebula as proxies to the Sun in order to 
determine standard abundances in the solar neighborhood? 

Table [5] lists the abundances of light elements in the Sun, Or ion B stars, and Orion neb- 



ula, showing an excellent agreement for all elements if we adopt iFeldman fc Widingj (120031 ) 
values for neon and argon in the solar photosphere. This table has been constructed from 
representative studies that we view as current state-of-the-art works in different objects, and 
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without intending to review comprehensively all abundance studies in the solar neighbor- 
hood (in particular for C, N, and O). We aim at putting the present study in a broader 
context and discuss its relevance to the Ar solar abundance question. First, we need to 
emphasize that the agreement between the nebular and the photospheric analyses for all five 
elements (C, N, 0, Ne, and Ar) cannot be fortuitous. It rather demonstrates the degree of 
reliability that can be achieved today in abundance analyses. We believe that we may thus 
reasonably claim that the accuracy of these abundances is 0.1 dex or better. Second, the 
long-standing puzzle regarding the C, N, O abundances in the S un and in young OB stars 



appears now to have been resolved. In lAnders Grevessd (119891 ) reference work, the CNO 



abundances were higher in the Sun than in OB stars, in an apparent contradiction with the 
olde r age of the Sun. Rec ent studies have now led to revise downward the solar abundances, 



and lAsplund et al.l (120051 ) recommend C, N, O solar abundances that are very close to the 
values derived for Orion stars and nebula. Therefore, these two arguments provide a robust 
support to the idea that the composition of Orion's young material (in stars and in the 
nebula) is very similar to the solar composition; hence, Orion's studies may be good proxies 
to determine the composition in the solar neighborhood, in particular for neon and argon 
too. 

Let us assume on the other hand that the coronal abundances are representative of the 
standard composition in the Solar system. This assumption thus implies that the chemical 
enrichment over the last 5 Gyr would have been a factor of 2 for neon, and a factor of 
2 to 3 for argon, while C, N, and O remained virtually unchanged. Analyses of PNe do 
not support such an enrichment of neon and argon over the last few billion years. Within 
1 kpc of the Sun, 7 PNe have mean abundances o f <A(0)> = 8.66, <A(Ne)> = 8.12 and 



<A(Ar)> = 6.50 (IPottasch fc Bernard-Salad 120061 ) . in good agreement with Orion's results 



(see Table [5]). The 7 PNe reveal however an enrichment in carbon and nitrogen resulting 
from nucleosynthesis in the late evolution of their parent stars. We therefore conclude that 
neon and argon do not reveal a strong chemical enrichment in the solar neighborhood during 
the last few billion years, and that the depletion of neon and argon in the solar corona most 
likely arise from yet little- understood processes in the corona, such as the FIP effect. 



5. CONCLUSIONS 

We have completed a new analysis of Ar II lines in the blue spectrum of 10 B main- 
sequence stars in the Orion association. The analysis involved a fully consistent NLTE 
treatment based on our NLTE line blanketed model atmospheres. We have derived a mean 
argon abundance which is in excellent agreement with recent results for the Orion nebula 
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( lEsteban et al.ll2004l ). However, the derived abundance is significantly higher (0.45 dex) tha n 
the standard solar value recommended in the recent compilation of lAsplund et al.l (120051 ). 
Based on a review of available argon abundance determinations in various sites in the solar 
neighborhood, we argue that the Orion's results much more likely represent the abundance 
of argon (and neon) in the solar neighborhood, while the lower coronal abundances reflect a 
depletion due to physical processes in the corona such as the FIP effect. 

Therefore, we propose to adopt the following values as new reference abundances in the 
neighborhood of the Sun in the Galactic disk: 



A(Ne) = 8.09 ± 0.06 
A(Ar) = 6.63 ± 0.10 



The argon abundance is in excellent agreement with Lodders (2008) recent work, in which she 
adopts for the solar system (protosolar) abundance, A(Ar) = 6.57 ± 0.10. These abundances 
yield the abundance ratios Ne/O = 0.26 and Ar/O = 0.009. We believe that these values 
should provide valuable reference points for various fields of investigations, from the solar 
corona to the ionization of the interstellar medium. 



We thank Martin Asplund for inquiring about the argon abundances in B stars, Carrie 
Trundle for sending equivalent widths of Ar II lines measured in the star HD 135485 and 
Letizia Stanghellini for discussions. The work reported here is supported in part by the 
National Science Foundation through AST03-07534, AST03-07532, AST06-46790 and NASA 
through NAG5-9213. 
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Table 1. Observing log. 



Star ID 


Spectral Type 


V 


UT Date 


Exposure times [s] 


HD 35039 


B2 IV-V 


4.70 


2007-01-30 


2x150 


HD 35299 


B1.5 V 


5.68 


2007-01-30 


2x300 


HD 35912 


B2 V 


6.38 


2007-02-06 


720 








2007-03-03 


2x700 


HD 36285 


B2 IV-V 


6.32 


2007-03-03 


2x600 


HD 36351 


B1.5 V 


5.46 


2007-02-06 


330 








2007-03-03 


500, 2x700 


HD 36591 


Bl IV 


5.33 


2007-02-06 


180 


HD 36959 


Bl V 


5.67 


2007-02-06 


2x270 


HD 37209 


Bl V 


5.72 


2007-02-06 


360 


HD 37356 


B2 IV-V 


6.18 


2007-03-03 


2x600 


HD 37744 


B1.5 V 


6.20 


2007-02-06 


480 
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Table 2. Selected Ar II lines and oscillator strengths. 



Wavelength [A] 


E cxc [cm x ] 




Accuracy 3, 


4426.001 


135 085.9960 


0.158 


A 


4430.189 


135 601.7336 


-0.174 


A 


4589.898 


148 620.1411 


0.100 


A 


4657.901 


138 243.6442 


-0.236 


B 


4726.868 


138 243.6442 


-0.103 


A 


4735.906 


134 241.7392 


-0.108 


A 


4764.865 


139 258.3384 


-0.06 


B 


4806.020 


134 241.7392 


0.210 


A 


4847.810 


135 085.9960 


-0.223 


B 


4879.864 


138 243.6442 


0.246 


A 


5062.037 


135 601.7336 


-0.465 


B 



a Estimated accuracy is better than 3% (A) or better 
than 10% (B). 
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Table 3. Argon and oxygen abundances. 



Star 




log 5 


A(Ar) 


A(O) 


HD 


35039 


20550 


3.74 


6.62 ± 0.06 


8.60 


HD 


35299 


24000 


4.25 


6.70 ± 0.04 


8.57 


HD 


35912 


19590 


4.20 


6.66 ± 0.05 


8.70 


HD 


36285 


21930 


4.40 


6.59 ± 0.04 


8.80 


HD 


36351 


21950 


4.15 


6.59 ± 0.05 


8.76 


HD 


36591 


26330 


4.20 


6.73 ± 0.09 


8.60 


HD 


36959 


24890 


4.40 


6.67 ± 0.07 


8.76 


HD 


37209 


24050 


4.15 


6.63 ± 0.04 


8.83 


HD 


37356 


22370 


4.15 


6.76 ± 0.04 


8.67 


HD 


37744 


24480 


4.40 


6.65 ± 0.05 


8.63 
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Table 4. Argon abundances in various sites. 



A(Ar) 



Sites, methods 



References 



6.66 ± 0.06 
6.50 ± 0.05 
6.54 ± 0.15 a 

6.62 ± 0.05 
6.6 - 6.7 

-0.14 



6.50 



-0.22 



6.1 ± 0.2 

6.62 ± 0.12 

6.45 ± 0.03 
6.23 ± 0.1 
6.38 



-0.18 
-0.30 



6.51 



-0.26 
-0.17 



6.18 ± 0.03 
6.26 ± 0.1 



B stars; NLTE 
B stars; NLTE 
B stars; LTE 

Orion nebula, t 2 = 0.022 
Orion nebula, t 2 = 0.0 - 0.35 

Planetary nebulae, d < 1 kpc 

ISM, Ar I AA1048, 1067 

Low-altitude impulsive flare, Ar/Mg 

Solar corona, X-ray flares, Ar XVII 
Solar corona, Ar XIV A4412 
Solar corona, X-ray flares 

PMS stars in Orion, corona, X-rays 

Solar energetic particles, Ar/O 
Lunar soils, Ar/O 



This study 

Holmgren et al. (1990) 
Keenan et al. (1990) 

Esteban et al. (2004) 

Peimbert fc Torres- Peimbert (1977) 

Pottasch fe Bernard-Salas (2006) 

Sofia & Jenkins (1998) : Lehner et al. f2003) 

Feldman & Widing (1990 , 2003) 

Phillips et al. (2003) 
Young et al. (1997) 
Veck fc Parkinson (1981) 

Maggio et al. (200Z) 

Reames 1998) 
Cerutti (1974) 



'Corrected as discussed in text. 



Table 5. Abundances in the solar neighborhood. 



Element Sun Orion B stars Orion nebula PNe (d < 1 kpc) References 



He 


10.98 ± 0.02 




10.988 ± 0.003 




2,..., 7, ... 


C 


8.39 ± 0.05 


8.35 ± 0.05 


8.42 ± 0.02 


8.65«; 2 ? 


1, 4, 7, 8 


N 


7.78 ± 0.06 


7.78 ± 0.07 


7.73 ± 0.09 


o i q+0.14 
°- io -0.22 


1, 4, 7, 8 





8.66 ± 0.05 


8.70 ± 0.09 


8.65 ± 0.03 


8.66l°;°? 


1, 4, 7, 8 


Ne 


8.08 ± 0.10 


8.11 ± 0.04 


8.16 ± 0.09 


o i O+0.14 
°- ±z -0.20 


3, 5, 7, 8 


Ar 


6.57 ± 0.12 


6.66 ± 0.06 


6.62 ± 0.05 


O.OU_ 22 


3, 6, 7, 8 



References. — (1) Asplund et al. ( 2005 ); (2) Christensen-Dalsgaardl ( 1998 ); 

(3) iFeldman fc WidineJ d2003h : (4 N ) ICunha fc Lambertl Jl994h: (5 ) Cunha et all kood \: (6) This 
paper: (7) festeban et all (bopj ): (8) IPottasch fc Bernard- Salasl (l2006h . 
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Fig. 1. — Ar IIAA4426, 4430 lines in the whole stellar sample. The stars in the panels are 
arranged in order of increasing T e g. 
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Fig. 2. — Sample Ar II lines in the spectrum of the B1.5V star, HD 35299: observed 
spectrum (blue points) and best-fit NLTE model spectrum (gray line; red line in electronic 
edition). The effect of changing the argon abundance by 0.1 dex on the predicted line profile 
is illustrated in the top left panel. 
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Fig. 3. — Top panels: NLTE ionization fractions of neon and argon as function of depth 
in a model atmosphere, T eff = 24, 000 K, log g = 4.25, and solar composition. Dashed lines 
show LTE ionization fractions for the same physical conditions. Bottom panels: Departure 
coefficients of Ne and Ar + levels of the analyzed lines (lower levels: full lines; upper levels: 
dashed lines); the thick lines show the departure coefficients of the Ne and Ar + ground states. 



